function activation = computeRBF(state, centre, width)
    activation = exp(-norm(state - centre)^2 / width^2);
end